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Abstract 

Using the the concept of cracking we explore the influence of density fluctuations and local anisotropy 
have on the stability of local and non-local anisotropic matter configurations in general relativity. This 
concept, conceived to describe the behaviour of a fluid distribution just after its departure from equi- 
librium, provides an alternative approach to consider the stability of selfgravitating compact objects. 
We show that potentially unstable regions within a configuration can be identified as a function of the 
difference of propagations of sound along tangential and radial directions. In fact, it is found that these 
regions could occur when, at particular point within the distribution, the tangential speed of sound is 
greater than the radial one. 

1 Introduction 

An increasing amount of theoretical evidence strongly suggests that a variety of very interesting physical 
phenomena may take place giving rise to local anisotropy, i.e. unequal radial and tangential stresses Pr Pi_ 
(see [1,2], and references therein). In the newtonian regime it has been pointed out in the classical paper 
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by J.H. Jeans [3], and in the context of General Relativity, it was early remarked by G. Lemaitre [4] that 
local anisotropy can relax the upper limits imposed on the maximum value of the surface gravitational 
potential. Since the pioneering work of R. Bowers and E. Liang [5] its influence in General Relativity has 
been extensively studied. 

Any model for an anisotropic compact object is worthless if it is unstable against fluctuations of its 
physical variables and, different degrees of stability/instability will lead to different patterns of evolution in 
the collapse of self-gravitating objects. Therefore, as expected, stability of anisotropic matter configurations 
in General Relativity has been considered since the beginning of the effort to understand the effects of 
tangential pressures on a selfgravitating matter configuration [5]. Very soon, in 1976, W. Hillebrandt and 
K.O. Steinmetz [6], considering the problem of stability of fully relativistic anisotropic neutron star models, 
showed (numerically) that there exists a stability criterion similar to the one obtained for isotropic models. 
Later, Chan, Herrera and Santos [7] studied the role played by the local anisotropy in the onset of dynamical 
instabilities. They found that small anisotropics might drastically change the evolution of the system. 
Recently, an analytical method has been reported to extend the traditional Chandrasekhar's variational 
formalism [8] to anisotropic spheres [9]. 

L. Herrera introduces, in 1992, the concept of cracking (or overturning) [10] which is a qualitatively 
different approach to identify potentially unstable anisotropic matter configurations. The idea is that fluid 
elements, at both sides of the cracking point, are accelerated with respect to each other. It was conceived 
to describe the behaviour of a fluid distribution just after its departure from equilibrium. Later on, Herrera 
and collaborators [7] showed that even small deviations from local isotropy may lead to drastic changes in 
the evolution of the system as compared with the purely locally isotropic case. More over, they found that 
perturbations of density alone, do not take the system out of equilibrium for anisotropic matter configura- 
tions. Only perturbations of both, density and local anisotropy induce such departures [11,12]. This concept 
refers only to the tendency of the configuration to split (or to compress) at a particular point within the 
distribution but not to collapse or to expand. The cracking, overturning, expansion or collapse, has to be 
established from the integration of the full set of Einstein equations. Nevertheless, it should be clear that 
the occurrence of these phenomena could drastically alter the subsequent evolution of the system. If within 
a particular configuration no cracking (or overturning) is to appear, we could identify it as potentially stable 
(not absolutely stable), because other types of perturbations could lead to its expansions or collapse. 

In the present paper we shall explore the influence that fluctuations of density and local anisotropy 
have on the possible cracking (or overturning) of local and non local anisotropic matter configurations in 
general relativity. We show that, for particular dependent perturbations, potentially unstable regions within 
anisotroic matter configurations could occur when the tangential speed of sound, dP±/dp, is greater than 
the radial, dPr/dp. This can give a more clear physical insight when considering the stability of particular 
anisotropic configurations when independent perturbations occurs. 

This paper is organized as follows. Section [2] will describe our notation through a brief discussion 
of local anisotropy matter configurations. The concept of cracking for selfgravitating anisotropic matter 
configurations and its relation with the sound speeds, is considered in Section [3) The models and modeling 
strategy are presented in sections |4] and [5l Finally some results and conclusions are displayed in Section [6l 

2 Anisotropic matter configuration in General Relativity 

We shall consider a static spherically symmetric anisotropic distribution of matter, described by the Schwarzchild 
line element ds^ = e^^'^^'dt^ — e'^'^^^dr^ — r^{d9^ + sm6d(f>^) and having an energy- momentum tensor repre- 
sented by = diag [p, ~Pr, —P±, —P±], where, p is the energy density, Pr the radial pressure and P± the 
tangential pressure. For this matter configurations, the general relativistic hydrostatic equilibrium equation 
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can be written [5] as 

dP,. , , /m + 47Tr^Pr\ 2 , ^ 

dr \ r (r — 2m) J r 

Obviously, in the isotropic case {P± — Pr) it becomes the usual Tolman-Oppenheimer-Volkov (TOV) 
equation, which constrains the internal equilibrium structure of general relativistic, isotropic, static perfect 
fluid spheres and it is considered in standard textbooks of gravitation [13, 14],. 

It is clear that the last term in[Tl (P^ — Pr) = represents a "force" due to the local anisotropy. 
This "force" is directed outward when P^ > P^ <^ A > and inward if Pj^ < P^ <^ A < 0. Therefore 
we should have more massive configurations if A > and less massive ones if A < 0. This becomes more 
evident when the most extreme situations i.e. Pj^ ^ and Pr = or Pj_ = and Pr ^ are considered. 



2.1 Ansatze for an anisotropic equation of state 

If a density profile, p — p{r) is given, it is possible to integrate [T] when the definition of mass and the two 
other equations of state, 

m(r) = 47r /" pf'^df, Pr = Pr(p) and P_i = P_i (P,.) , (2) 
Jo 

are provided. The first equation of state, P,- — Pr{p), corresponds to the standard barotropic equation of 
state for time-independent systems. In order to close the system, we have also to provide a second equation 
of state relating radial and tangential pressures. Pi = P± (Pr)- It has been shown [15,16], that there exists 
a unique global solution to ([T]) if: p is a continuous positive function; Pj_(r) is a continuous differentiable 
function; Pr{r) is a solution to the equation with starting value P_l(0) = Pr(0), and both pressures are 
positive at the center (i.e. P±(0) = Pr(0) > 0), therefore in the whole interior of the body. 

Much of the efforts to disentangle the physics of very dense matter is reflected by the various "radial" 
equations of state: P^ = Pr{p) availables (see [17,18] and references therein). In contrast, very little is known 
for the much less intuitive second equation of state Pj_ = P± (Pr)- This is the reason why different ansatze 
are found to introduce anisotropy in matter configurations (see, for instance, references [1,2,5,19-28]). The 
unknown physics in the "tangential" equation of state is partially compensated by using heuristic criteria 
(geometric, simplicity or any other assumption relating radial and tangential pressures). Therefore, most of 
the exact solutions for the differential equation ([T]) found in the literature have been obtained from excessively 
simplifying heuristic assumptions and, in addition, some of the conditions to become "physically acceptable 
fluids" are not verified. 



2.2 Acceptability conditions for anisotropic matter 

The interior solution should satisfy some general physical requirements. Some of the "physical acceptability 
conditions" for anisotropic matter have been stated elsewhere [1,2] as 

1. density, p, radial pressure, Pr, and tangential pressure. Pi, should be positive everywhere inside the 
configuration. 

2. gradients for density and radial pressure should be negative, 

dp „ . dPr „ 

/ < , and -^<0- 
or or 

3. inside the static configuration the speed of sound should be less than the speed of light, 

— — < 1 and — — < 1 ; 
op op 
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4. in addition to the above intuitive physical requirements, the interior solution should satisfy [29] either: 

• the Strong Energy Condition: p + P,. + 2 P± > 0, p + Pr > and p + P± > or 

• the Dominant Energy Condition: p > Pr and p > P± 

5. junction conditions [30], match the matter configuration to the exterior Schwarzchild solution. Because 
of the continuity of the First Fundamental Form, the definition of mass in[2l evaluated at the boundary, 
becomes the total mass, M ~ m{a), as measured by its external gravitational field. More over, the 
continuity of the Second Fundamental Form forces the radial pressure to vanish at the boundary, r = a, 
of the sphere Pr\r=a — ^■ 

Notice that as a consequence of the junction conditions, the radial pressure should vanish at the boundary, 
but not the tangential one. However, both should be equal at the centre of the matter configuration. Also 
notice that there is no restriction on the gradient for the tangential pressure. 

These reasonable physical requirements validate the assumptions made for both equations of state and, in 
many cases, exclude possible mathematical solutions of the system [T] and [5] Delgaty and Lake [31] considering 
the isotropic case {Pr — P±), found that, from 127 published solutions only 16 satisfy the above conditions. 
In particular, it is worth mentioning that, in order to have a causal theory of matter we have to demand 
that the sound speed be, at most, the speed of light. This important requirement when obviated, violates 
physical principles usually required for a matter physical theory [32]. 



3 Instability and cracking of anisotropic compact objects 

As we have stressed above, in a series of papers Herrera and collaborators [10-12] elaborated the concept 
of cracking for selfgravitating isotropic and anisotropic matter configurations. It was introduced to describe 
the behaviour of fluid distributions just after its departure from equilibrium, when total non- vanishing radial 
forces of different signs appear within the system. 

This section will describe the general framework of the cracking approach to identify potentially stable 
(and unstable) anisotropic matter configurations. We explicitly use some of the tacit assumptions for mod- 
eling cracking (or overturning) within these matter configurations and, finally, we propose a more intuitive 
criterion based on the difference of sound speeds to estimate the relative magnitude for the density and 
anisotropy perturbations and to evaluate the stability of bounded distributions. 



3.1 Cracking: the general framework 

Herrera and collaborators state that there is cracking whenever the radial force is directed inward in the 
inner part of the sphere and reverses its sign beyond some value of the radial coordinate; or, when the 
force is directed outward in the inner part and changes sign in the outer part, we shall say that there is an 
overturning. These effects are related to the tidal accelerations of fluid elements [11,33], defined by 



ds ds 



(3) 



where 5x'^ is a vector connecting the two neighbouring particles; h'^ denotes the projector onto the three-space 
orthogonal to the four- velocity and du"/ds = u^^u"^. More over, defining 



dPr , , f ra + inr^Pr 

ar \ r [r — Ira) 



r 



(4) 
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it can be shown that [3] and IH evaluated immediately after perturbation, lead to [11,12] 



R = - 



e\p + Pr) 



W2~2de 

ar e ' r — — , 
ds 



(5) 



where & represents the expansion. Again, here, ds^ — e^'^'^^dt^ — e'^^^'-'dr^ — r^{d9^ + sinddcfP), the static 
Schwarzchild line element, has been assumed (see reference [11] for details). [4] is just the hydrostatic equi- 
librium equation [1] that vanishes for static (or slowly evolving) configurations. It can be appreciated from 
([5]) that for cracking to occur at some value of < r < a, it is necessary that dQ/ds vanishes somewhere 
within the configuration. It is also clear the non local nature of this effect; and that small deviations from 
local isotropy may lead to drastic changes in the evolution of the system as compared with the purely locally 
isotropic case [12]. 

3.2 Cracking revisited 

Following [12], we assume that the system having some pressure and density distributions satisfying R — 0, 
is perturbed from its hydrostatic equilibrium. Thus, fluctuations in density and anisotropy induce total 
radial forces (i? ^ 0) which, depending on their spatial distribution, may lead to the cracking; i.e radial force 
directed inward, (i? > 0), or, overturning, directed outward, (R < 0) of the source. Therefore, we will be 
looking for a change of the sign of R, beyond some value of the radial coordinate. We will exclusively consider 
perturbations on both, density and local anisotropy, under which the system will be dynamically unstable. 
In other words, Sp and SA are going to be considered as independent perturbations; but fluctuations in mass 
and radial distribution pressure depend on density perturbations, i.e. 



'""^-Sp, 



p + Sp 



Pr{p + Sp,r)^Pr{p,r)+SPr^Prip,r) , ^ 
m{p + Sp, r) = 47r ^^{p + 6p)f^df k, m{p, r) + ^r^6p 



Now on, expanding [4] we have, formally. 



r. r. , r. AN ^R, OR OR^ OR^^ 

R « Ro{p, Pr, m, A, r + —dp + tt^SP,- + i^Sm + —5A 

op oPr dm oA 



(6) 



(7) 



and by using [6] it can be shown that 

R = dp 

where 



dR 



47r ^dR 
3 dm 



2SA 

r Sp 



dR 
dp 



m 



in Prr^ 



> and 



dR _ {p + Pr) {1 + 8tt Py) 



d' 



'm 



{r -2 my 



> 



(8) 



(9) 



r(r — 2m) 

It is immediatelly seen that, in order to have R = and consequently a change in its sign: 

• both, the anisotropy and the density, have to be perturbed; 

• both, anisotropy and density perturbations, have to have the same sign, i.e. SA/Sp > 0. 

In order words, potentially stable configurations should have SA/Sp < everywhere because R never changes 
its sign [34]. 



5 



3.3 Cracking and sound speeds 

When arbitrary and independent density and anisotropy perturbations are considered (as in all previous 
works concerning cracking [10-12,34]) there is little physical criteria to establish the size (absolute and/or 
relative) of the perturbation, i.e. how small (or big) the perturbations should be. Different orders of 
magnitude (and relative size SA/Sp) of perturbations could produce a cracking but we could be describing 
an unphysical scenario. Additionally, all these previous works only consider constant perturbations. It is 
possible that variable perturbations could be more efficient inducing cracking within a particular matter 
configuration. Again, there is no criteria in establishing the functionality of the perturbation throughout the 
matter distribution. 

We are going to consider a particular type of dependent perturbation whose relative order of magnitude 
could be bounded by the behavior of some physical variables and could be checked by physical intuition. Ob- 
viously, general perturbations should be independent because they emerge from non related physical phe- 
nomena. But we are looking for some physical variables whose behavior could be checked in order to identify 
potential cracking. 

It is easy to convince oneself that 



SA S (Pl - Pr) SPj_ SPr 2 2 

1^ Ji^""'^ '''''' ' 



where v^^. and v'^j_ represent the radial and tangential sound speeds, respectively. 

This will be the key concept, we will use in revisiting Herrera's approach to identify potentially unstable 
anisotropic matter configurations based on the concept of cracking. Now, by considering the sound speeds 
and evaluating I10|, we could not only have a more precise idea of the relative order of magnitude of the 
perturbations {SA and Sp) but also what are the regions more likely to be potentially unstable within a 
matter configuration. 

It is clear that, because < < 1 and < v'^j_ < 1, we have jw^j^ — v'^rl — 1- Thus, 

-1 < "J^sX — '>^sr Potentially stable , 

1 < V^s± -vl<l^{ (11) 

< v'^^_ — < 1 Potentially unstable . 

Therefore, we can now evaluate potentially unstable regions within anisotropic models based on the differ- 
ence of the propagation of sound within the matter configuration. Those regions where u^^ > v'^j_ will be 
potentially unstable. On the other hand, if v"^^ < everywhere within a matter distribution, no cracking 
will occur. It is worth mentioning, concerning this criterion, one of the extreme matter configurations men- 
tioned above {P± ^ and Py. = 0) is allways potentially stable for cracking; and the other one (Pi = 
and Pj. 7^ 0) becomes potentially unstable. 

More over, for physically reasonable models, the magnitude of perturbations in anisotropy should always 
be smaller than those in density, i.e. \v'\^ — ^srl — ^ ^ l*^^! — \^P\- When 5A/5p > 0, these perturbations 
lead to potentially unstable models. 

Next section will be devoted to explore the effectiveness of this criterion on the stability of bounded matter 
configurations, having different equations of state. Again, we recall that the concept of cracking refers only 
to the tendency of the configuration to split and its occurrence has to be established from the integration of 
the full set of Einstein equations. In addition, it is clear that there could also be some perturbations that 
do not induce cracking but could cause instabilities that lead the configuration to collapse or to expand. 
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4 Perturbations and cracking for anisotropic configurations 



In order to illustrate the workability of the above criterion we shall work out several density profiles 
for models satisfying the physical acceptable conditions. Thus, in addition to the positivity of density and 
pressures profiles, their gradients and the fulfillment of the energy conditions (strong or dominant), we shall 
pay special and particular attention to the conditions bounding sound speeds (radial and tangential) within 
the matter configuration. 

The idea will be to provide the density profile; then to obtain the radial pressure, Prif) from a "radial" 
equation of state Pr — Pr{p{r)) and next to solve the tangential pressure P±{r) from the anisotripic TOV 
m That is. 

r dPr (p + Pr) fm + inr^Pr 



Pir)'^Pr = PMr))^P.-Pr + -^^ + ^[-^^;-^) ■ (12) 

Then, the radial and tangential sound speeds are calculated and their difference \v'^±^ — Vg^\ is evaluated. 
Next, by using [TT] the pontential stability or unstability is established. This will be confirmed by a change 
in the sign of R described by [5] and [5] 

To illustrate the above criterion [11] we shall analyze four cases concerning qualitatively different density 
profiles. We have selected two local (one singular and one non singular) and two non-local conformally flat 
anisotropic solutions. By local models we mean the standard way to express an equation of state where the 
energy density and radial pressure are related at a particular point within the configuration, i.e. P = P{p{r)). 
On the other hand, by non local models we will understand those where the radial pressure Pr (r) is not only 
a function of the energy density, p{r), at that point; but also its functional throughout the rest of the 
configuration. Any change in the radial pressure takes into account the effects of the variations of the energy 
density within the entire volume [27,35,36]. It has been shown that in the static limit, this particular radial 
equation of state can be written as 

Pr{r) = p{r) - ^ [ f2p(f)df. (13) 
Jo 

It is clear that in equation (fT3l) a collective behavior of the physical variables p{r) and Pr{r) is present. 
4.1 Anisotropic Tolman VI model 

This model was introduced by Cosenza, Hcrrera, Esculpi and Witten [20] starting from the singular Tolman 
VI density profile [37]. The original isotropic Tolman VI solution is not deprived of a physical meaning. 
It resembles a highly relativistic Fermi Gas with the corresponding adiabatic exponent of 4/3. By using a 
heuristic method these authors determine other physical variables representing an anisotropic static matter 
configuration; i.e. 

K 3 /l-^/^^ 3 /21-25,/^\ 



where, the junction conditions force the adjustment of the parameter, to -R" = 3/567r, and the radius is given 
by a = 81/49. 

Sound speeds can be determined from 1141 and can be written as 

7(7 + 3^-9./^) , 3 (49 + 25^^-70./^) 

(7 -3 VI) 4(7-37^) 
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By using [T3] the particular expresion for [5] can be obtained for this model as 



Rto. 



ImanVI 



2Sp 
r 



(588+ 180^- 6727^) 



vi ^ + 539 + 195f 



658 



16(7-37^)^ 



SA 



(16) 



It is worth mentioning that this model does not fulfill all the acceptability conditions stated in 12.21 It is 
singular at the center and near the boundary surface (r / a ~ 0.706) the tangential pressure becomes negative. 
Despite this unphysical situation, this model is presented because, as it will later become clear in Section 
O the difference in sound speed is constant through the whole configuration, which represents the above 
mentioned constant perturbation relation considered in previous works [10-12,34]. 

4.2 Non local Stewart Model 1 

This model emerges from a density profile proposed by B. W. Stewart [21], to describe anisotropic conformally 
flat static bounded configurations; which was also, recently, considered for non local anisotropic matter 
distributions [27]. Starting from this density profile we can find Pr{r) and P±{r) as 



P = 



Pr - 



1 (e^^*^ - 1) (e*^'^ + mre'^^'- - 1) 



87rr2 



1 (l-e-^-^'-)(e4^'--8AW^''-l) 
87rr2 



Pi = 



7r[l + e2-'^'-]4 



(17) 
(18) 
(19) 



Again, the parameter K has to be obtained from the junction conditions; which means that K has to 
satisfy a trascendental equation 



1 



g4Ka_g^^g^«a_i^Q ^ K=—hi 

2a 



From [TTI [18] and [19] the corresponding sound speeds can be found 



(20) 







"(e2^'-2)'-3 




(g4Kr _ 1)2 _^ 8^2^2g2ifr 


(g2Xr _ 2)^ - 3 





and 



32if3^3g4Kr (g2Kr _ i) 


(e2^'- + l)"' 


(g4Kr _ 1)2 ^ 8^2^2g2/fr- 


(g2ifr _ 2)^ - 3 



Now, equation [8] can also be obtained for this model as 

25p 



R 



N LStewartl 



[{e^^-" + 1) (5 + vl^) + AKr {e^^-^ - l)] Kr SA 
6(e2^'- + l)^(e2^'--l)"^ 'Sp 



(21) 



(22) 



(23) 
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4.3 Non local Stewart Model 2 

This is a second density profile proposed by B. W. Stewart [21] which was recently proved to be non local [36] 



1 



Pr = - 



87rr2 
1 

87rr2 



1 - 



sin(2ifr) sin^(ii:r) 



Kr 



Pi - 



Kr 
1 

8^ 



1 - 



sm"^ (Kr) 



(24) 
(25) 
(26) 



As in the previous models, junction conditions (M ~ m{a) and Pr{a) — 0) determine the coupling 
constant K . In this case it has to satisfy also a trascendental equation 



svaKa 
Ka 



2M 



1/2 



a I 



and cos Ka 



1 



3M 



1 _ 2M 



Thus, 



K 



Fromf 



a~2M 

] and [55] the corresponding sound speeds can be found 



fi\n^ [Kr) (3 - K"^ 



I sin(2Xr)i^r 

and 



Kr 
2 



Kr 



i(2Jfr) 



siT?{Kr) 



[cofi{Kr)Kr — sin(iirr)] 
Now, equation [8] can also be obtained for this model as 



[cos{Kr)Kr — sin(_ftrr 



Rnls 



tewart2 



2Sp 



2 sin(2i^r)i^r - 9 sin^ (i^r) (w^^ + l) SA 
12 sin^(i^r) [sin(2i^r)i^r - 2 sm^{Kr)] 



(27) 



(28) 



(29) 



(30) 



4.4 Florides-Stewart-Gokhroo & Mehra Model 

This density profile is due originally to P.S. Florides [38], but also corresponds to different solutions, con- 
sidered by Stewart [21] and, more recently, by M. K. Gokhroo and A. L. Mehra [23]. The Florides-Stewart- 
Gokhroo-Mehra solution represents densities and pressures which, under particular circumstances [39] , give 
rise to an equation of state similar to the Bethe-Borner-Sato newtonian equation of state for nuclear mat- 
ter [33,40,41]. 
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P = P = Pc[l 



Kr^ 



5-3K 



5-3K 



15/ie' 



-2A 



j (5 - 3K) 



5w 
5-3K 



1 - -Kt]-' ] (1 - Kt]-' 



with 



M _^ 2fi7j^ (5 - 3Ktj^) r 

^=— , e =1 — — '- and ?? = -. 

a 5 — a 



and pc the density at the center of the matter configuration. 

From [Hll [32] and [33] the corresponding radial sound speed can be found as 

2 _ 2^ (5 - QKrf) [1 - (1 + n) if] +n[b-iK{l + 2^7;" 



Kj (5 - 2,K) (1 - 7^2)^ 



and 



2 2 



QpKrf [2-rf{2 + n) 
(5 - m) (1 - 7^2) 



1 — n o — A 



'7^ \ , »7~" 



2e 



-A 



where 



30/xe-^ (l - rf 



j(5-3i^) 



2n-qe-^ \ 4:njKf] (4 - 3X772) 
1 — 772 



1 



5-3K 



and 



5-3/C 

15^i 2Tr{5-3K) 
~r~ {1- 772)"+! 

4^77 (5 - 6/^772) 



47re 



-A 



7^6 



-A 1 



1-772 

f,J (5 - 3j^7/2) 
(1 -i^772)-^ 



5-3K 

For this example equation [5] can also be obtained as 



R 



FSGM 



25p 



2^1-2 



01^, ■ -2A 15e-^ (1 - 77^) - J (5 - 6/^772 

(5 — 3iv)je 2Aq^ \ ' \ 

2/x (5 - 3^^772) [l5e-^ (1 - 772)" + j (10 - 9^^772)] 
f + 25) + 30j [e-^ (1 -Tf-)K- j] 77 



6/z77^ 



(5-3ifrj2 



-75e-2Ml-,72)'" 



(5A 
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Density Profile 

TolmanVI 
NL Stewart 1 
NL Stewart 2 
Gokhroo & Mehra 



M/a 
0.21 
0.32 
0.39 
0.26 



M{Mq) 
1.42 
2.15 
2.68 
1.76 



■a 



Pa X 



IQ^^ {gr / CTTT') xlO^^ {gr/cm^) 



0.31 
0.65 
1.19 
0.44 



2.30 NA 

6.80 1.91 

8.49 2.14 

0.00 2.09 



Table 1: All parameters have been chosen to represent a possible compact object with a = 10 Km. and the 
corresponding mass function satisfying the physical acceptability and energy conditions 

5 The modeling performed 

As it can be appreciated from the parameters displayed in table [1] all the models considered have radius 
(a = 10 Km.) and total masses, M (in terms of solar mass Mq) that correspond to typical values for 
expected astrophysical compact objects. The boundary redshifts Za, surface and central densities, pa and Pc 
that emerge from our selection, also fit the typical values for these objects. 

Profiles for the radial, w^^, and tangential, sound speeds, as well as its difference, v'^j_ — w^^, are 
displayed in Figure [TJ The perturbation relation, SA/Sp = v"^^ — v^^ fulfills the physical restriction — 1 < 
SA/6p < 1 for all models considered. Notice that SA/6p, is constant within the matter distribution for 
the Tolman VI anisotropic model (plate / in Figure [1]) . This type of constant perturbation relations were 
standard for modeling cracking in previous works [10-12,34]. Because SA/Sp < 0, the sound speed stability 
criterion, [TTl states that in the Tolman VI anisotropic model, no cracking will occur. Non local Stewart 
models are sketched in plates // and ///, respectively. For these two models we could implement variable 
perturbation relations, SA/Sp, through the matter configuration; because —1 < SA/Sp < 0, no cracking 
will occur in these models either. Finally, the most interesting scenario emerges from the Florides-Stewart- 
Gokhroo-Mehra model [23] with j — 7, K — 1 and n = 2, shown in plate IV. As it is evident from this plate, 
the perturbation relation, SA/Sp not only has a variable profile, but it also changes its sign, alternating 
potentially stable and unstable regions within the distribution. In fact, this model presents two potentially 
unstable regions: < = r/a < 0.2570 and 0.7565 < = r/a < 1 where SA/Sp > 0. 

The profiles of R for each model are ploted in Figure [2] and the above stability assumptions can be 
contrasted with the change in sign for the expressions \TE[ [23l [30] and [39l It is clear from the i?-plots 
displayed in this figure, that the models of Tolman VI, NL Stewart 1 and NL Stewart 2 models do not 
present any cracking point (plates /, // and ///, respectively). On the other hand, the Florides-Stewart- 
Gokhroo-Mehra model displays a cracking point at 77 w 0.17986 within the first potentially unstable region 
,0<V = r/a< 0.2570. 

6 Results and conclusions 

We have revisited the concept of cracking for selfgravitating anisotropic matter configurations introduced by 
L. Herrera and collaborators [10-12]. It has been shown that for some particular dependent perturbations, 
the ratio for fluctuations in anisotropy to energy density, SA/Sp can be interpreted in terms of the difference 
of sound speeds, i.e. SA/Sp ~ — u^^; where Vg^ and v'^j_ represent the radial and tangential sound 

speeds, respectively. It is evident from [TT] that regions where > v-^j_ will be potentially unstable. On 
the other hand, if < w^j^ everywhere within a matter distribution, no cracking will occur and it could be 
considered as stable. 

This reinterpretation could be useful to refine and make the concept of cracking more physically related 
to the potential instability due to the behavior of some physical variables within matter configurations. It 
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Figure 1: Variations of the radial and tangential sound speeds for anisotropic configurations. 
Plates I, II, III y IV represent Tolman VI, NL Stewart 1, NL Stewart 2 and Gokhroo & Mehra, 
respectively. Curves a, b and c correspond to v^^., j_ y j_ — v"^^, respectively. 
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Figure 2: The induced total radial force R for the anisotropic configurations. Plots /, //, /// 
and IV represent Tolman VI, NL Stewart 1, NL Stewart 2 and Gokhroo & Mehra, respectively. 
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is easy to determine each sound speed, their difference [TU] and the sign of the difference. Thereafter, we 
could clearly identify from [8] which regions are more likely to be potentially unstable within a particular 
matter distribution. This can be appreciated from the Florides- Stewart- Gokhroo-Mehra model (Figure 
Opiate IV) which displays a cracking point at 77 « 0.17986 within the first potentially unstable region 
0<T] = r/a< 0.2570. 

Additionally, because each sound speed has to be less than the speed of light, it implies that their 
difference has the physical restriction: |(5A/(5p| ^ \'^'s± ~ ^sr\ — 1- This is very important in order 
to characterize a particular model as potentially unstable. It is possible to find cracking points within a 
configuration for unphysical set of fluctuations in anisotropy and energy density, i.e. |(5A/(5p| > 1, but the 
existence of these cracking points could not lead to physical potentially unstable models. More over, the 
physical restriction, \5A/Sp\ < 1 also conditions the relative order of magnitude of the perturbations. 

The ratio of perturbations, SA/Sp are now not necessarily constant. Models considered in previous 
works [10-12, 34] have constant fluctuations, because there were no other criteria to introduce in order 
to evaluate the change in the sign of R. Now, the possibility to introduce variable fluctuations based on 
difference of sound speeds, enrich the applicability of the cracking framework to evaluate instabilities within 
anisotropic matter configurations. 

It is worth mentioning that, concerning this criterion, one of the extreme matter configurations mentioned 
above Pi 7^ and Pr = is always potentially stable, and the other P± = and Pr could experiment 
a cracking (or overturning) scenario. The study of matter configurations with vanishing radial stresses traces 
back to G. Lemaitre [4] and for non static models have been considered in [1]. Recently, this model has 
been studied [42,43] concerning its relation with naked singularities and conformally flat models has been 
considered in [24]. Extreme models with vanishing tangential stresses seems to be useful describing highly 
compact astrophysical objects having very large magnetic fields {B > 10^^ G) [44] 

As we have pointed out, any model for a static compact object is worthless if it is unstable against 
fluctuations of its physical variables. If a particular static model is unstable against these fluctuations it could 
follow different possible patterns in its subsequent evolution. It could collapse, expand, split or overturn. 
Perturbations play a crucial role not only evaluating the stability of a particular static model, but identifying 
trends in possible future evolution of the model. Their study should be considered from different points of 
view and formalisms. In this work we have considered only those perturbations, related through radial and 
tangential sound speeds, that lead to identify potentially unstable regions. Independent perturbations (not 
related via any physical quantity) could also exist and could also lead to cracking (or overturning) points but, 
in this case there is no criteria to quantify their order of magnitude. Other types of perturbations leading to 
expanding or collapsing evolutions could be considered in the standard Chandrasekhar's variational formalism 
( see [8,9] and references therein). Again, we stress the fact that those different possible evolution patterns for 
unstable configurations, refers only to a tendency. Its occurrence has to be established from the integration 
of the full set of Einstein equations. 
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